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We prove that the electron density function of a real physical system can be uniquely determined 
by its values on any finite subsystem. This establishes the existence of a rigorous density-functional 
theory for any open electronic system. By introducing a new density functional for dissipative 
interactions between the reduced system and its environment, we subsequently develop a time- 
dependent density-functional theory which depends in principle only on the electron density of the 
reduced system. 
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Density-functional theory (DFT) has been widely used 
as a research tool in condensed matter physics, chemistry, 
materials science, and nanoscience. The Hohenberg- 
Kohn theorem [jj lays the foundation of DFT. The 
Kohn-Sham formalism provides a practical solution 
to calculate the ground state properties of electronic sys- 
tems. Runge and Gross extended DFT further to calcu- 
late the time-dependent properties and hence the excited 
state properties of any electronic systems 0- The accu- 
racy of DFT or time-dependent DFT (TDDFT) is de- 
termined by the exchange-correlation functional. If the 
exact exchange-correlation functional were known, the 
Kohn-Sham formalism would have provided the exact 
ground state properties, and the Runge-Gross extension, 
TDDFT, would have yielded the exact time-dependent 
and excited states properties. Despite their wide range 
of applications, DFT and TDDFT have been mostly lim- 
ited to isolated systems. 

Many systems of current research interest are open 
systems. A molecular electronic device is one such sys- 
tem. DFT-based simulations have been carried out on 
such devices 0, IE IE S IE 0, E3, E3, E3 These sim- 
ulations focus on steady-state currents under bias volt- 
ages. Two types of approaches have been adopted. One is 
the Lippmann-Schwinger formalism by Lang and cowork- 
ers 0. The other is the first-principles noneq uilibrium 
Green's function (NEGF) technique g |E EE 111 El • In 
both approaches the Kohn-Sham Fock operator is taken 
as the effective single-electron model Hamiltonian, and 
the transmission coefficients are calculated within the 
noninteracting electron model. The investigated systems 
are not in their ground states, and applying ground state 
DFT formalism for such systems is only an approxima- 
tion E3- DFT formalisms adapted for current-carrying 
systems have also been proposed recently, such as Kosov's 
Kohn-Sham equations with direct current .14] and Burke 
et a/.'s Kohn-Sham master equation including dissipation 
to phonons E3- However, practical implementation of 
these formalisms requires the electron density function 
of the entire system. In this paper, we present a rigor- 
ous DFT formalism for any open electronic system. The 
first-principles formalism depends in principle only on the 
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electron density function of the reduced system, and can 
be used to simulate the transient electronic response. 

It is often implicitly assumed that the electron density 
function of any real physical system is real analytic |lfi| . 
This is reflected in many existing first-principles method- 
ologies. In practical quantum mechanical simulations, 
real analytic functions such as Gaussian functions, Slater 
functions and plane wave functions are adopted as basis 
sets, which results in real analytic electron density dis- 
tribution. However, analyticity is not guaranteed for all 
presently used basis functions. For instance, in the lin- 
earized augmented plane wave (LAPW) method E3 > ra ~ 
dial functions times spherical harmonics span the inside 
of muffin-tin spheres whereas plane waves are used for 
the interstitial regime. By this approximation the elec- 
tron density cannot be analytically continued across the 
interface of the two disconnected parts. It is thus in- 
teresting to ask whether the electron density function of 
any real system is in principle real analytic. This ques- 
tion was settled recently for time-independent systems. 

Fournais et al. proved that the electron densities of 
atomic and molecular eigenfunctions are real analytic in 
r-space away from the nuclei |18| . Their proof is based on 
the fact that the time-independent Schrodinger equation, 
= , is an elliptic partial differential equation, and 
its eigenfunctions can always be real analytic (except iso- 
lated points) and quadratically integrable. Their rigor- 
ous proof implies that the electron density function of 
a real physical system is real analytic (except at nuclei) 
when the system in its ground state, any of its excited 
eigenstates, or any state which is a linear combination of 
finite number of its eigenstates. In their derivation [18j, 
nuclei are treated as point charges, and this leads to non- 
analytic electron density at these isolated points. Note 
that the isolated points at nuclei can be excluded for 
the moment from the physical space that we consider, 
so long as the space is connected. Later we will come 
back to these isolated points, and show that their in- 
clusion does not alter our results. We show now that 
for a time-independent real physical system the electron 
density distribution function in a sub-space determines 
uniquely its values on the entire physical space. This 
is nothing but analytic continuation of a real analytic 
function. The proof for univariable real analytical func- 
tions can be found in textbooks, for instance, Ref. E3- 
The extension to multivariable real analytical functions 
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is straightforward. 

Lemma : The electron density function p(r) is real 
analytic in a connected physical space U. D C U is a 
sub-space. If p(r) is known for all r £ D, p(r) can be 
uniquely determined on the entire space U. 

Proof: To facilitate our discussion, the following nota- 
tions are introduced. Set Z + = {0, 1,2,.. .}, and 7 is an 
element of (Z+) 3 , i.e., 7 = (71,72,73) G (Z+) 3 . The dis- 
placement vector r is denoted by the three-dimensional 
variable x — (xi,X2,x^) G U. Denote that 7! = 
71 ! 72 !73 ! , x-y = xf x? xf, and £^ = ^ g? 2 Q 3 . 

Suppose that another density distribution function 
p'{x) is real analytic in U and equal to p{x) for all x G D. 

We have ^ = for all x G D and 7 G (Z+) 3 . 

Taking a point x Q £ D and at or infinitely close to 
the boundary of D, we may expand p{x) and p(x') 



around Xq, 
and p'(x) = V 
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that the convergence radii for the Taylor expansions of 
p{x) and p'(x) at xo are both larger than a positive fi- 
nite real number 6, we have thus p{x) = p'(x) for all x G 



d~<p(x) 
dxi 



d-'p'jx) 
dxi 



Db(xo) — {x : \x — xq\ < b} since 

Therefore, the equality p'{x) — p{x) has been expanded 
beyond D to include -D&(xo). Since U is connected the 
above procedures can be repeated until p'(x) = p{x) for 
all x 6 U . We have thus proven that p can be uniquely 
determined in U once it is known in D. 

With the above Lemma we are ready to prove the fol- 
lowing theorem: 

Theorem 1: A connected time-independent real phys- 
ical system is in the ground state. The electron density 
function p(r) of any finite subsystem determines uniquely 
all electronic properties of the entire system. 

Proof: Assuming the physical spaces spanned by the 
subsystem and the connected real physical system are D 
and U, respectively. D is thus a sub-space of U, i.e., 
D C U. According to the above lemma, p{r) in D deter- 
mines uniquely its values in U , i.e., p(v) of the subsystem 
determines uniquely p(r) of the entire system. 

Inclusion of isolated points, lines or planes where p(r) 
is non-analytic into the connected physical space does not 
violate the above conclusion, so long as p(r) is continuous 
at these points, lines or planes. This can be shown clearly 
by performing analytical continuation of p(r) infinitesi- 
mally close to them. Therefore, p(r) of any finite subsys- 
tem determines uniquely p{v) of the entire physical sys- 
tem including nuclear sites. Hohenberg-Kohn theorem Q 
states that the ground state electron density distribution 
of any system determines uniquely all its electronic prop- 
erties. Hence we conclude that ground state p(r) of any 
finite subsystem determines all the electronic properties 
of the entire real physical system. 

As for time-dependent systems, the issue is less clear. 
Although it seems intuitive that the electron density 
function of any time-dependent real physical system is 



real analytical (except isolated points in space-time), it 
turns out quite difficult to prove the analyticity rigor- 
ously. Fortunately we are able to establish a one-to- 
one correspondence between the electron density function 
of any finite subsystem and the external potential field 
which is real analytical in both i-space and r-space, and 
thus circumvent the difficulty concerning the analyticity 
of time-dependent electron density function. For time- 
dependent real physical systems, we have the following 
theorem: 

Theorem 2: The electron density function of a real 
physical system at to, p(r,to), is real analytic in r-space, 
and the corresponding wave function is $(^0)- The sys- 
tem is subjected to a real analytic (in both i-space and r- 
space) external potential field v(r, t). Let D be a finite r- 
subspace. The time-dependent electron density function 
on the subspace D, p{r,t) with r G D, has thus a one-to- 
one correspondence with v(r,t) and determines uniquely 
all electronic properties of the entire time-dependent sys- 
tem. 

Proof: Let v(r, t) and v'(r, t) be two real analytical po- 
tentials in both t-space and r-space which differ by more 
than a constant at any time t ^ to, and their correspond- 
ing electron density functions are p(r,t) and p'{r,t), re- 
spectively. Therefore, there exists a minimal nonnegative 
integer k such that the fc-th order derivative differentiates 
these two potentials at to- 
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[v(r,t)-v'{r,t)} 



7^ const. 
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Following exactly the Eqs. (3)-(6) of Ref. we have 

Qk+2 



g t k+: 



[p(r,t)-p'(r,t)} 



-V-«(r), (2) 



t=*o 



where 



u(r) = p(r,i )V< 



dt k 



[v(r,t)-v'{r,t)] 



(3) 



Due to the analyticity of p(r,to), v(r,t) and v'(r,t), V • 
u(r) is also real analytic in r-space. It has been proven 
in Ref. that it is impossible to have V ■ u(r) — on 
the entire r-space. Therefore it is also impossible that 

V • u(r) — everywhere in D (otherwise based on the 
Lemma proven earlier, V-u(r) can always be analytically 
continued from the particular finite r-subspace where its 
values are to the entire r-space, and this would lead to 

V • u(r) = everywhere in the entire r-space). We have 
thus 



Qk+2 

dt k + 2 



[p{v,t)~ p'(T,t)] 



* 



(4) 



for r G D. This confirms the existence of a one-to-one 
correspondence between v(r, t) and p(v,t) with r G D. 
p(r, t) on the subspace D thus determines uniquely all 
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FIG. 1: Schematic representation of the experimental setup 
for quantum transport through a molecular device. 



electronic properties of the entire system. This completes 
the proof of Theorem 2. 

Note that if <5(io) is the ground state, any excited 
eigenstate, or any state as a linear combination of finite 
number of eigenstates of a time-independent Hamilto- 
nian, the prerequisite condition in Theorem 2 that the 
electron density function p(r, to) be real analytic is auto- 
matically satisfied, as proven in Ref. 18]. As long as the 
electron density function at t = to, p(r,to), is real an- 
alytic, it is guaranteed that p(r, t) on the subsystem D 
determines all physical properties of the entire system at 
any time t if the external potential v(r, t) is real analytic. 

According to Theorem 1 and 2, the electron den- 
sity function of any subsystem determines all the elec- 
tronic properties of the entire time-independent or time- 
dependent physical system. This proves in principle the 
existence of a rigorous DFT-type formalism for open elec- 
tronic systems. All one needs to know is the electron den- 
sity of the reduced system. The challenge that remains 
is to develop a practical first-principles formalism. 

Fig. n depicts an open electronic system. Region D 
containing a molecular device is the reduced system of 
our interests, and the electrodes L and R are the envi- 
ronment. Altogether D, L and R form the entire system. 
Taking Fig.^as an example, we develop a practical DFT 
formalism for the open systems. Within the TDDFT for- 
malism, a closed equation of motion (EOM) has been de- 
rived for the reduced single-electron density matrix o~(t) 
of the entire system po): 



ia(t) = [h(t),o-(t)], 



(5) 



where h(t) is the Kohn-Sham Fock matrix of the entire 
system, and the square bracket on the right-hand side 
(RHS) denotes a commutator. The matrix element of a 
is defined as o~ij(t) = (a.j(t) ai(t)}, where a,i(t) and ctj(i) 
are the annihilation and creation operators for atomic 
orbitals i and j at time t, respectively. Fourier trans- 
formed into frequency domain while considering linear 
response only, Eq. © leads to the conventional Casida's 
equation |2l|. Expanded in the atomic orbital basis set, 
the matrix a can be partitioned as: 



o~l old o- lr 
o~dl o~d &DR 

0~RL 0-rd &R 



(0) 



where <tl, an and <jd represent the diagonal blocks cor- 
responding to the left lead L, the right lead R and the de- 
vice region D, respectively; cfld is the off-diagonal block 
between L and D; and ctrd, <tlr% &dl, o~dr and grl 
are similarly defined. The Kohn-Sham Fock matrix h 
can be partitioned in the same way with a replaced by h 
in Eq. @. Thus, the EOM for cjd can be written as 



i&D = [hD,o-n]+ ^2 {hoaPaD ~ o- Da h aD ) 

a=L,R 

= [h D ,(T D ]-i ^ Qa, (?) 



a=L,R 



where Ql (Qr) is the dissipative term due to L (R). 
With the reduced system D and the leads L/R spanned 
respectively by atomic orbitals {1} and single-electron 
states {k a }, Eq. (7J is equivalent to: 



icr,. 



leD 



a=L,R 



AS) 

(9) 



k a £a 



where m and n correspond to the atomic orbitals in re- 
gion D; k a corresponds to an electronic state in the elec- 
trode a (a = L or R). h n k a is the coupling matrix el- 
ement between the atomic orbital n and the electronic 
state k a . The current through the interfaces Sl or Sr 
(see Fig. ^| can be evaluated as follows, 



Ja(t) 



"// r ^ (r ' t} 

= - E jt^v 

k a £a 

— *E E ( fl k a i°~ik a - o-k a i hk a ) 

= -5^Qa,B = -tr[Q a (*)], (10) 

leD 

i.e., the trace of Q a . 

At first glance Eq. JHJ is not self-closed since the dissi- 
pative terms Q a remain unsolved. According to Theorem 
1 and 2, all physical quantities are explicit or implicit 
functionals of the electron density of the reduced system 
D, po{r,t). Note that pi)(r, t) = p(r,t) for reD. Q a 
is thus also a universal functional of pu(r,t). Therefore, 
Eq. JHJ can be recast into a formally closed form, 



ia D = 



h D [r,t; p D (r,t)],a D 



a=L,R 



) a [r,t; p D (r,t)}. 



(11) 

Neglecting the second term on the RHS of Eq. 
leads to the conventional TDDFT formulation in terms 
of reduced single-electron density matrix [2(} for the 
isolated reduced system. The second term describes 
the dissipative processes between D and L or R. Be- 
sides the exchange-correlation functional, an additional 
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universal density functional, the dissipation functional 
Q a [r,t; p D (r,t)], is introduced to account for the dissi- 
pative interaction between the reduced system and its 
environment. Eq. {TTJ is the TDDFT EOM for open elec- 
tronic systems. Burke et al. extended TDDFT to include 
electronic systems interacting with phonon baths 
they proved the existence of a one-to-one correspondence 
between v(r,t) and p(r,t) under the condition that the 
dissipative interactions (denoted by a superoperator C in 
Ref. 0|) between electrons and phonons are fixed. In 
our case since the electrons can move in and out the re- 
duced system, the number of the electrons in the reduced 
system is not conserved. In addition, the dissipative in- 
teractions can be determined in principle by the electron 
density of the reduced system. We do not need to stip- 
ulate that the dissipative interactions with the environ- 
ment are fixed as Burke et al.. And the only information 
we need is the electron density of the reduced system. 
In the frozen DFT approach |22( an additional exchange- 
correlation functional term was introduced to account for 
the exchange-correlation interaction between the system 
and the environment. This additional term is included 
in hjj[r,t; pi)(r,t)] of Eq. An explicit form of the 

dissipation functional Q a is required for practical imple- 
mentation of Eq. lfTT|l . Admittedly Q Q [r, t; /3d (r, t)] is an 
extremely complex functional and difficult to evaluate. 
As various approximated expressions have been adopted 
for the DFT exchange-correlation functional in practical 
implementations, the same strategy can be applied to the 
dissipation functional Q a . Work along this direction will 
be published elsewhere [23| . 

Given Q a [p] how do we solve Eq. in practice? 

Again take the molecular device shown in Fig. ^ as an 
example. We may integrate Eq. (|llfl directly by satis- 
fying the boundary conditions at Sl and Sr. The only 
boundary condition we need is the potentials at Sl and 



Sr. We need thus integrate Eq. together with a 
Poisson equation for Coulomb potential. And the Pois- 
son equation is subjected to the boundary condition de- 
termined by the potentials at Sl and Sr. It is impor- 
tant to point out that although in principle its physical 
span can be small, in practice the reduced system is to 
be chosen so that Eq. can be solved readily with 
convenient boundary conditions. For instance, for the 
molecular electronic device depicted in Fig. ^ the re- 
duced system D contains not only the molecular device 
itself, but also portions of the left and right electrodes. 
In this way the Coulomb potential at the boundary take 
approximately the values of the bulk leads. 

To summarize, we have proved rigorously the existence 
of a first-principles method for both time-independent 
and time-dependent open electronic systems, and devel- 
oped a formally closed TDDFT formalism by introduc- 
ing a new dissipation functional. This new functional Q a 
depends only on the electron density function of the re- 
duced system. With an explicit form for the universal 
dissipation functional Q a , the time evolution of an open 
electron system in external fields is fully characterized 
by Eq. (|llfl . In practical calculations, we need thus fo- 
cus only on the reduced system with appropriate bound- 
ary conditions. This work greatly extends the realm of 
density-functional theory. 
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